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ENGLir:i ABSTRACT 

The formation of primordial black holes (PBH) in the case 
of an expanding Universe was discussed under the assumption 
that at the early stages Uie expansion is near to quasi-isotropic. 
The choice of a quasi-isotropic solution is brought about by the 
fact that, even though quasi-isotropy does not exist, the particle 
generation near the singularity should inevitably result in a 
quasi-isotropic solution. 

For simplicity the spherically symmetrical case will be con- 
sidered below. Then the quasi-isotropic solution is character- 
ized by one function only. 

The problem is solved by means of a computer. The computa- 
tions performed allow the following conclusions to be made: 

1. PBH form when there are important deviations from the 
Fridmanian model. If the deviation from the Frldmanian model 
can be imaged as a part of the closed world that is smoothly 
matched with the flat world, the degree of the deviation from 
the Frldmanian model can be quantitatively characterized by the 
Lagrange radius of the mentioned matching. If the Lagrange 
radius is assumed to be the unit at which the maximum of the 
Euler radius is achieved, our computations show that for the 
black hole formation it is necessary that the radius of the 
matching will be 0.85 to C.9. 

2. The hydrodynamlcal computations made the fact clear that 
the smoothness of the matching is of great Importance; this can 

be quantitatively characterized by a relative width of the match- 
ing A with reference to the Lagrangian radius The smaller A, 

the greater the pressure gradients preventing contraction. 
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Ouyr computaMo^ns show ^that J>he„ pjt’esauce^^j^^ prevents - - 

PBH formation decreasing the mass of the hole as compared with 
that, which should have occurred under the same initial condi- 
tions but with zero pressure. 

3. To plot the total hydrodynamical picture of accretion 
it is necessary to compute several dozens of hydrodynamical 
times. It will be done in another paper. But it is seen even 
here that the similarity solution typical for catastrophic ac- 
cretion cannot be obtained. It is explained by the matter 
ejection from the disturbed region due to pressure gradients; 
as a result, the size of a black hole, if it can nevertheless 
form, will be much smaller than the horl zon and the formation 
moment . 


HnM)DMMlGS OP PRIMORDtAL BLACK 

D.K. Nadezhln, ‘I.D. Novikov, A/G. Polnarev 

1. Introduction /5* 

In 1966> Ya. B. Zel’dovlch and I. D. Novikov [l-3l> and 
then in 1971 Hawking [if] pointed out the possibility of the for-^ 
mation of black holes at the very beginning of cosmic expansion 
from primary matter. Later on, msny works discussed the problem 
of primordial black holes (PBH) (see, for example, [5-7]). 

Interest in PBH Increased particularly after Hawking's dis- 
covery [8] of quantum evaporation of small mass black holes be- 
cause PBH can have small masses. The Hawking process is Important, 
on the one hand, for the physics of the early stages of expansion 
of the Universe and, on the other hand, as a possible method of 
detecting PBH in today's Universe [9-12]. 

The following two problems are basic for the theory of PBH: 

1) lifhat must the deviations be from the Fridman cosmic 
model at the beginning of expansion so that PBH will form? 

2) How will accretion of matter on the PBH formed occur? 

Both problems were posed and considered in the very first 
works on PBH [1,2,3]. In these and succeeding works, it was 
pointed out that the complete answer to both questions required 
a numerical calculation on a computer. Besides PBH in literature, 
white holes are considered as well. The Importance of quantum 
processes close to singularity and the process of accretion [1-3] 
converting a white hole into a singular black hole was pointed 
out [ 13 ] for these objects. /6 


^Numbers in the margin indicate pagination in the foreign text. 


This article clarifies a description of the method and re- 
sults of calculation of the hydrodynamics of the phenomena indi- 
cated using an EVM in the simplest hypothesis as to the spherical 
symmetry of the processes considered. 

2. Posing the Problem 

PBH cannot form in the early stage of expansion of the Uni- 
verse (with the equation of state P = e/3)» if deviation from the 
Fridman model was small both in density and in metrics. Actually, 
as Lifshits has pointed out [l4], perturbations of density which 
are small in ampll^de 6e/c increase only to the values of pertu- 
bation oi^ metrics which arh assumed to be small. After this, the 
condition of increase of perturbatrdRB'hi^ea^? 

scale of perturbations becomes smaller than the 'horCajgrr/^;^^^^ they 
return to sonic oscillation. 

On the other hand, as is well known, deviations from the 
Fridman model exist which result in the formation of PBH. Actually, 
let us consider the nearest singular state of a uniform semiclosed 
world [ 15 - 17 ], matched through a narrow orifice with a uniTorm flat 
cosmic model. During evolution of a semiclosed world, the signal 
from the orifice moving with the speed of sound successfully 

reaches the internal sections of the semiclosed world where these 
parts will already be at the s_tage of contraction under an intrinsic 

gravitational radius and consequently, as is well-known-, they form 
a black hole. 

The question arises: what must the critical value of de- 

viation from the Fridman model be during which a black hole forms 
and what is the smallest deviation where the hole does not form. 

This problem is solved in our work. It will be formulated more 
precisely later on. We are also interested in the hydrodynamic 
processes which accompany the formation of PBH. /7 
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Thus, we will consider some deviations from the Fridman model 
iLn the ini t tar state . SeTectidn wh^^^ Initial conditions 

of close singularity very large in a general case (see [18,19]). 
However, it is important to note that the quantum processes of 
close singularity sharply limit the possible initial conditions. 

As has been pointed out in [20,21], immediately after completion 
of quantum phenomena (that is, after the moment t = tp^ = VhG/c® 
the solution of the OTO*- equation must be quasl-isotropic [22,23]. 
This solution is characterized by the fact that close to slngu- 
,larity, all of the elements of matter expand according to the same 

'on 

law e = 32 n Gt^ (for P = ^ e), as is true in the Fridman model. 

At the same time, the spatial cross section t = const can be a 
random distortion method. Assignment of this distorted cross 
section is a general Initial condition for a quasi-isotropic 
solution. In the case we considered of spherical symmetry of a 
curve df .-three dimensional space, it is described by a single 
functlpn, of the radial coordinate. Thus, we will consider the 
cohdiJt'ion for PBH occurring depending on assignment of this func- 
tion. 


The process of formation of PBH will be considerec? depend- 
ing on- the amplitude of deviation from the Fridman model of close 
singularity (origin of expansion) and to the profile of this ^ 
deviation (however, staying within the framework of spherical 
symmetry). We will prescribe a deviation close to singularity 
in the form of a spherical field with the accompanying space 
of a constant positive curve — that is, the perturbed field 
corresponds to part of a closed Fridman model. The amplitude 

^General theory of relativity. 

2The time of formation of a black hole with mass M and the time /]_ 
during which this black Jiole then vaporizes, involve an order 
of magnitude of the relationship ty^p/tpi,,^ = where 

- 5*10"^d — the Planck time. From this it is apparent that 
^vap 'tfripj lO'^g (that is, the Planck mass). Con- 

sequently, in hydrodynamic calculation, we correctly do not 
consider the quantum processes of vaporization of black holes. 
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of deviation can be described with the aumber whlch^ sbaracterl ses 
the part of the closed space of a constant positive curve which 
we have taken (see paragraph 4). This part of a closed world 
through a transition field is matched to the flat Fridman model. 

The width of the transition field, as our calculations have shown, 
will be the second significant parameter of the problem. 

Here It is assumed that within the perturbed field the solu- 
tion is precisely described by the Fridman model. Thus, pertur- 
bation is such that the full mass of substance within the perturbed 
field is precisely equal to the mass which would occur in this 
field in an unperturbed Fridman model. 


3 . Initial Equations 

As has been pointed out above, we will consider a spherical 
problem. We will write the spherically symmetric metric 

= C^e^cLi^ - cL(f^ ) . ( D 


Let us select the accompanying system of calculation. Then the R- 
Lagrange radius. Here we will describe the Einstein system of 
equations for the metric (1) (see [24]). We note only that this 
system has a first interval found in the work of Podurets [25] 


and along with the laws of conservation (T^.j^ 
the following simpler system: 


= 0) equivalent to 




m 


I 


HT 


(3) 


(4) 





Here m — is mass Included within the Lagrange radius R, e = dens- 
ity of energy, P is pressure. 

Numerical calculation of the system of equations (2-6) as 
Podurets proposed [25] appears to be Ineffective for the problems 
which we have selected to solve. The fact is that during forma- 
tion of black holes in an expanding Universe, in the cent'-’al field 
of perturbation, expansion of matter is -changed by compression. 
Near to this moment in time (different for different R) r = 0 and 
equation (2) has a root characteristic. Inasmuch as this moment 
was unknown earlier, it appears more convenient to increase the 
magnitude of the system (2-6) and obtain a new system free of 
root characteristics. 


Let us convert the system (2-6). Incidentally, we will in- 
troduce dimensionless values of r, m, e, connected to r, t, 
m and e in the following way 



Here r is the gravitational radius corresponding to character- / 

s 

istlc mass. Later on we will omit the symbol remembering 
that we are working with dimensionless values. 


Let us Introduce the value u related to f by the relation- 


u -- has a ^hyaicai valoslty oT expansion or contrac- 
tion of the Lagrange surface R = const. 

Prom 2-6), differentiating (2) in time, we obtain the fol- 
lowing system: 









u 


0 




(9) 

( 10 ) 


£ 
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( 11 ) 


Equation (8) (see also [26]) is a relativistic analogy of 
the second law of Newton. On the left is acceleration and on 
the right — force. The first member is force involving the 
gradient of pressure and the two succeeding members describe 
gravitational attraction. 

We will assume that pressure P depends only on e (this 
hypothesis is not true for a model of stars but is completely 
suitable for the early stages of expansion of the Universe), and 
then the equation of state can be written In the form 


( 12 ) 


(13) 


Integrating (5) and (6), we have; 


Here 






^ =e»p[| 



r(f)£ 


/II 


(14) 


(15) 


If Y does not depend on e, then A = f(t) and j>(R) — are 

random functions of time and radius. 


We will consider the formation of PBH in a cosmic Fridman 
model with a flat accompanying space. Let us select; the function 
f(t) in such a way that our accompanying system will have large 
values of R (where the solution is Fridmanian) simultaneously 
synchronous , that is 




( 16 ) 


It is not difficult to point out that for a uniform spatially 
flat Fridman solution with random values of y (not depending for 
simplicity on e) we have: 



with calculation of (7). ( 17 ) 


Then from (13) and (15-17) we find f(t): 

= — )C^/rhM) (18 

Selection of the function (j»(R) will be considered in the 
next paragraph. 
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Later on we will have to determine the portion of matter 
entering into a black hole. For this it is necessary to find 
at each moment of time the boundary of the black hole — the 
horizon of tne event. As is well known (see [6]), a local 
criterion for the position of the horizon of the event does not 
exist and to find it, it is necessary to integrate zero geodesics 
in the future. In our problem, where the fall of matter into 
the black hole is important, the position of the horizon of the 
event at this moment essentially depends on the matter which will 
fall into the black hole in the future. For an observer who is 
near a black hole, it is more significant to determine the bound- 
ary of the black hole as the visibility horizon, that is, to de- 
termine the Lagrange radius corresponding to mass which at a 
given moment is under its gravitational radius. In this work we 
will use this determination of the boundary of a black hole. This 
boundary can be defined from calculation of local properties. 
Actually (see for example [6]), this boundary corresponds to the 
condition ^ = 1, and moreover, the requirement that the accompany- 
ing system be compressed into this boundary. 

Thus, we will form the final criterion of a black hole thus 
(using the dimensionless values defined above). 


If, for a certain point with a Lagrange radius R, two condi- 
tions are simultaneously fulfilled 

m ^ i 

( 19 ) 


'i 


> i 

LL < j 


then this point lies under the horizon cf visibility of the black 
hole. The equation of this horizon of visibility is 



( 19 ' ) 


8 


We note that the horizon of the event always lies Jtflthin the 
Tjorizon of vicibiHty and, generally speaking, in direct proxlm 
Ity to it. 

Initial and Boundary Conditions 


As was discussed in paragraph 2 we consider that near sin- 
gularity, expansion occurs in the quasi-lsotr^pic way. In the 
case of spherical symmetry this means that the soluvion has only 
one random function (ji(R) (see (1^)). An analytical solution of 
the system (8-11, 12, 13, 1^) or an equivalent system (2-4, 12- 
14) when t -► 0 (in singularity) with quasi-lsotroplc conditions 
gives us the following asymptotic solution 



f J J 



Here two members of expansion are considered. Near singularity 
the spherically symmetrical deviation from the Fridman model with 
a flat accompanying space occurs in the field 0 < R < R_^. Th;: 
center of deviation coincides with the origin of the coordinate. 


Using expansion (20-21), the Ir ‘ial conditions are pre- 
scribed in the following way; 


When R < R 


/] 4 
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( 22 ) 



UnH 


This corresponds to a space of a constant positive curve (that 
Is, a closed Fridman model). For a closed world, R changes from 
0 to It. Then the maximum of #(R) Is reached when 


When R > R^ = R_ + A, 

W‘ : 

I ), 


(23) 


this corresponds to a spatially flat Fridman solution. C — js 
a certain constant determined by the conditions of matching (22) 
and (23) in an intermediate field R_ < R < R^. 

When matching (22) and (23) we require continuity of the 
third derivatives of function <>(R), Inasmuch as these derivatives 
enter into equation (8) (in P* e', see (21) and in r*). 


As was pointed out in paragraph 2, we will select perturba- 
tion in such a way that outside the field of pertubatior. when 
R > R^, the solution would be precl.'?ely Prldmanlan, Therefore, 
it is obvious that in the transltioii field of matching R_ < R < R^ 
^ there must be a certain discharge (field of decreased density) 

so that the excess of mass within -vill be compensated in com- 
parison with a spatially flat Fridman solution. 


We will select time t^ fairly small, as the moment of time 
for prescribing the initial conditions so that the discarded 
members of expansion in (20) and in (21) will be small in com- 
parison with the remaining. Let us evaluate the value of the 
remaining ana eliminated members. By selection constant C can 
lead to smoothness of function ^(R), that is. 
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Then the maximum value of the second derivative In the order of 
magnitude equals 

.4 

— 7~~~~ 


(fy 


max 


(25) 


The first discarded member includes the third derivative produc- 
ing a value on the order of magnitude equal to 


ip** 

Y - * 


( 26 ) 


Prom this we find the condition for selection of the initial 
moment of time depending on the width of matching A: 




(27) 


Thus, with the consideration of everything presented above, 
the initial conditions, and this means the entire solution, are 
characterized by two values: 


1 ) 


6 = «-//? 


max 


3lR^ 


j , 


the ratio of the radius of matching to half the radius of a 
closed Fridman; 

2) A — is the width of matching. 

B characterizes the value of perturbation of geometry. The 
role of the pressure gradient depends on the value of A. With 
small values of A the role of the pressure gradient increases. 


The fact that the point is fairly distant (large values of 
R) we will use as the first boundary condition; the wave of re- 
lief did not successfully reach this and expansion occurs ac- 
cording to the law (17). The boundary condition on the left 
boundary also has a very simple form: 
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(28) 


Boundary conditions (28) would remain true at the latest 
stages of collapse after formation of a black hole when at the 
center for a finite time an actual physical singularity occurs. 
The moment of occurrence of singularity is different for dif- 
ferent Legrange radii, that is, singularity is not quasi-uniform 
which it was at first. Prom equations (2-6) or (8-11) one can 
obtain an asymptotic solution (Kazner type) which in the accom- 
panying system of computation with our selection of random func- 


tions has the form: 






(29) 


(30) 


xirhere 


n r 


f-3#i 




3 (S-f) 


’ f > Va 


(32) 


Thus, after singularity has formed, the left boundary con- 
dition continues to be at a certain Lagrange radius R, different 
from zero. This radius R is determined from the condition so 
that numerical solution is close to an asymptotic solution (29- 
31). During further collapse, as the new elements of matter drop 
to singularity, the position where the left boundary condition re- 
mains moves toward a large Lagrange radius. 


12 




■ 5. A .Method of Numerical Calculation 


For riumeriGal solution of a system of equations (8-11) the 
difference method was used similar to the basic characteristics 
in the methods used earlier for calculating gravitational col- 
lapse [26-27]. Therefore, we will not present here; the entire 
system of different equations but will limit ourselves only to 
a brief description of its basic peculiarities. The entire in- 
tegral of change of the Lagrange coordinates 0 x R < was 
broken up into 150 integrals and then the step of the break up 
&R was made smaller in the internal part and particularly close 
to the matching location of a closed world with the Fridman and 
was smoothly increased in the external part of the integral of 
change of R considered. In the difference equations, a method 
was used of "whole” and "semlwhole" calculation units according 
to space and time [28]. In the whole units, according to space, 
u, r, m were calculated and in the semiwhole — P,e and depending 
on them and A; moreover, u was calculated in semiwhole units 
according to time, and all of the remaining values — in whole 
units according to time. 



Each alternating step in time began with a calculation of u 
from the difference analog of equation (8), then from equation 
(9) and (10) sequentially distributions were calculated accoi ilng 
to space of r and m, and finally, the density of energy was found 
from the difference equation: 









which corresponds to (13); here J — is the number of the spatial 
calculation interval, and n — is the number of the time calcula- 
tion interval; J1 = ol50 — the whole number of space calculated 
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intervals. The values of density of energy In the flatter semi- /l8 
whole is defined -from. Ihe J^ovmdary w^ 

The difference system considered here was tested in sample 
calculations of gravitational collapse close to free fall, in 
calculations or Fridman expansion and Fridman expansion with 
small perturhations in the form of sonic waves. Then it was clear 
that the process of gravitational collapse is reproduced very 
well by a numerical method, that is, the different calculations 
of perturbation rapidiy attenuate as the central density increases. 
We note here than in the work of May and White [26] and Schwarts 
[273 only gravitational collapse was studied numerically. However, 
the series of calculations of the Fridman expansion (including the 
perturbed) showed that this numerical reproduction requires a 
more accurate approach inasmuch as the absolute value of calcu- 
lated perturbations of density decrease somewhat more slowly 
than the density itself and therefore relative errors gradually 
increase. So that one can calculate the stage of expansion 
fairly far, it is desirable to significantly decrease the step 
of integration in time in comparison with its maximum allowable 
value determined by the condition of numerical stability of 
Kurant . 


We used the following requirement as the criterion of 
selection of a variable step of integration according to time 
At: the ratio Ae/e does not exceed the prescribed fixed value 

C << 1 in all calculated spatial integrals (Ae — is change of 
e for one step in time). If, after completion of the alternating 
step of integration, A e/e > y appeared, although for one spatial 
interval, then At decreased and the subsequent step was recalcu- 
lated. It is convenient to use the Integral of equations (8-11) /19 
for control of precision: 


A . ~cj • i 

u 't “ 



(33) 


(see equation (2)). The difference of this expression from zero 
by several percentage points in relation to the component making 
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_the niaximurn contribution to the JLeft part (33) can be consi4ered^ 
fully exemplary. It was apparent that for the class of problems 
considered in this work this condition is fulfilled when ? < 2 * 

10”^. The results presented here of calculations were obtained 
using the values ^ = 10“'* - 10“^. 

6. A Description of Results of Calculation 

The results of calculation of the process of formation of ^ 

PBH depend on two parameters — amplitude of deviation from a 
flat Fridman model characterized by R_, and the width of matching 

A. 

As has already been discussed, in the field of matching 
there is a certain rarefaction. In the opposite case, outside 
there would not be a flat Fridman solution and one would have a 
perturbation of the ^ type where AM — is the excess mass (see 
in more detail [6]). Therefore, it is obvious that relative per- 
turbation of density when R_ < R < R_j_ is larger than the matched. 

As calculation showed, in the case of a very narrow matching, 
sometimes shock waves occur and solution of the problem in the 
end requires the introduction of viscosity. Figure 1 shows the 
evolution of this perturbation in the case where the width of 
matching A = 0.1R_, R_ = 1. Later on we will concentrate our 
attention on those cases where shock waves do not occur. Let us 
first select a certain fixed ratio of the width of matching to 
the radius of matching; let us say ^ =0.5 and we will study 

the following question: at what values of R_ does the black hole 
form and at v/hat values does it not form. It is clear that the 
smaller the values of R_, the smaller the perturbation metric is /20 
and with a fairly small R_, the black hole, as is well known, 
cannot form and perturbation must be transformed to sonic oscil- 
lation . 
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In Figure 5, R_ = that is, it amounts to 75? of 


Xt^ i^^apparent fi*om the Ctrawing tWt^^^ process occurs in 
the following way. At first, a general cosmic expansion occurs. 
The internal field expands more slowly that the external. Prom 
the field of matching to the center of perturbation and outside 
it there are waves of rarefaction. At a certain moment in time, 
in the internal fi 'Id, a shift in expansion to contraction occurs. 
During contraction, the pressure gradient increases so much that 
the central core is ejected and a wave of compression occurs out- 
side. In this case the perturbation of the metric, or, in other 
words, the gravitational field, was fairly large and the formation 
of a black hole did not take place. Perturbation was changed to 
sonic oscillations. 

In Figure 3 also a case is shown where a black hole does 
not form. Here R_ = perturbation of the metric is 

larger in this case and therefore the large and maximally allow- 
able density of energy occurs in the center of perturbation. The 
ratio of density of energy at the center of perturbation to dens- 
ity of energy outside perturbation reaches a maximum of 10^, that 
is, the problem is essentially nonlinear. Later on, perturbation 
is converted more and more to sonic waves. 


i } 


I 


F 


A qualitatively different picture is shown in Figure ^ 

(R = 0.85R_„„). We see here unlimited compression of the core. 

“ luaX 

The central densit / of energy tends toward the infinite. In 
order to be sure in this case a black hole actually forms, let 
us turn to Figure 5. In this drawing the ratio ^ is shown. In 
those fields where expansion occurs, r is increased and m, due 


to the forces of pressure, decreases. 


Therefore, the value p 


drops with expansion. When expansion changes to compression. 


then obviously the ratio p increases. 


If the pressure gradient 721 


m 


does not successfully become compression at the moment w'hen p 
reaches 1 (dashed line uji Figure then a black hole unavoidably 
occurs inasmuch as, in this case, both of the conditions fur the 
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criteria of a bla^ck hole are present (see paragraph 3 > formula 
(19)). The moment of formation of a black hole in our dimen- 
sionless units equals t = 2.39. The ratio of the mass of the 
black hole at the moment of formation to the mass of the entire 
perturbed field at this moment equals 20%. 

Thus, we can conclude that when a/R_ = 0.5> the black hole 

is formed only in the case where mass of a perturbed field is 

comparable to half the mass of a closed world, that is, R_ is 

close to R„„„ = The quantitative result of this is 
max 2 


when 

In Figures 6 and 7, similar relationships are presented 
R 

for a case when = 0.9. At this time, gravitational forces 

max 

even with a certain reserve exceed the forces of the pressure 
gradient. Therefore, the formation of a black hole occurs much 
earlier: t = I.87. And the ratio of the mass of a black hole 

to the mass of a perturbed field is larger at this moment than in 
the preceding case and equals approximately 30%. The latter in- 
dicates the fact that the force of pressure increased from the 
central field with a smaller mass than the preceding case. 

Now we will go on to study the question of how significant 

the matching width is for formation of a black hole. Let us 

take p -- - = 0.85, that is, the minimum R_ at which a black hole 
max 

forms and a decreased width in matching to A = 0.3R_. Then the 
role of the pressure gradient Increases and as is shown in Figure 
8, the black holes do not form. The pressure gradient success- /22 
fully leads to dissipation of the core earlier than a black hole 
is successfully formed. 




( 3 ^) 


= 0.5. 
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However, when R_ = °-9Rniax ^ " 0-3R_ and A = 0. 2R 


_ > 


the 


of 


black hole Is formed as when A = 0. 5R_> but the moment of forma- 
tion of the black hole ensues later than when A » 0.5R.J the ratio 

the mass of the hole formed to the mass of the perturbed rep;ion is smaller 
than in the cases indicated above. (When ' = 0.9, the black 

' max 

hole does not form even when A = O.IR ,) 


All of the calculated variations are presented together in 
Figure 9. The radius of matching is applied along the axis of 
the abclssa and the width of matching along the axis of the ordin- 
ate; the corresponds to the formation of a black hole, and ”0” 
to the occurrence of sonic oscillation. The dashed curve indi- 
cates the approximate boundary of the field of the parameters 
indicated where black holes occur. 


The Table presents the most important quantitative results 
for certain variations of calculation when a black hole formed. 


Thus, the calculations made make it possible to draw the 
following conclusions. 


PBH form only with very large deviations in the Fridman 

model corresponding to « 0.85-0.9. The width of the transl- 

max 

tion A affects the formation of PBH. The smaller A, the more 
significant is the role of pressure gradients, making it difficult 
for PBH to form. 


Before making numerical calculations, the role of pressure 
in forming PBH should be evaluated by constructing stationary 
[3] or selfmodeling [6] solutions. The hypothesis was presented 
that pressure can facilitate accretion of gas by the by the 
occurrence of PBH and result in a significant Increase in their 
mass. Carr and Hawking [6] pointed out that a selfmodeling 
solution does not exist which leads to catastrophic accretion 




- t;' 
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tjf the PBfr“inatter-dttr±ng-»irirt^3h“iftmeii5lon5 thre PBH^ tircreaBe /2H 
like the cosmic horizon. Our calculations Indicate that pres- 
sure actually interferes strongly with the formation of PBHj 
decrfeasing the mass of PBH in comparison with the mass which 
existed in the PBH under the initial perturbation but with the 
absence of pressure P =0. Actually, on the spatial section 
t = const near singularity, density of energy in the perturbed 
field R < R_ is larger than away from the Fridman model and 
the gradient of pressure force on R directed outward attempts 
to eject the mass. In the transition field R < R < R^, 
density is minimum and on its external boundary R^ pressure gradi- 
ent is directed Inward causing accretion. However, this phenome- 
non is less significant during the formation of PBH than the 
gradients described above for R_ which leads to "efflux" of the 
mass from the field of perturbation. For constructing a complete 
hydrodynamic picture of accretion on PBH (without being limited 
by the hypothesis of selfmodeling as in [6]), it is necessary to 
calculate the process fairly far in time on the order of lOrg/s. 

This we will do in a separate work. But right now one can say 
during the formation of PBH its dimension is much smaller than 
the cosmic horizon (the mass of the black hole at the moment of 
formation amounts to a total of only a few percentage points of 
the mass included within the cosmic horizon (see Table). There- 
fore, as was pointed out already in [1,2], catastrophic accretion 
does not occur and this conclusion coincides with the analysis of 
[5], made for selfmodeling solutions. 

The authors wish to thank Ya. B. Zel'dovlch for his interest 
in our work and his critical comment and also V. N. Lukash for 
his many comments . 
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TABLE 




PHYSICAL CHARACTERISTICS OP THE PROCE^SS-QE J^ORMATIOH OP 





0,5 

0 

1 

2 
3 

0,01 

0,62 

2539 

3,17 

2,31 

0,^ 

0,33 

0,42 

6,13 

1,01 

0,77 

0,74 

0,02 

1,12 

4,91 

6,30 

1 

0,00 ! 

0,00 
0,16 

0,32 i 


0 

0,01 

Z,42 

6,45 

0,02 

0,00 1 

0,5 

I 

0,62 

0,51 

1,09 

1,16 

0,00 .. 




2 

1,86 

0,48 

0,88 

3,79 

0,26 


3 

2,76 

0,49 

0,84 

5,55 

0,46 


0 

0,01 

2,42 

4,74 

0,02 

0,00 

0,3 

I 

0,62 

0,52 

0,86 

1,16 

0,00 


2 

2,42 

0,45 

0,67 

30 

0,14 


3 

3,12 

0,44 

0,65 

■-,29 

0,28 


0 

0,01 

2,42 

4,05 

0,02 

0,00 

0,2 

I 

0,62 

0,54 

0,78 

1,16 

0,00 







2 

3,31 

0,42 

0,56 

6,59 

0,09 


3 

3,80 

0,41 

0,55 

8,26 

0,15 


In the tabulated material are equivalent to decimals. J 


— mass within the Lagrange radius R_ (mass of part of a 
closed Prldman world); 

— mass within R^ (mass of total perturbation) ; 

— mass Included within the cosmic horizon; 

— mass of a black hole (mass Included within the visibility 
horizon) ; 

— initial moment; 1=1 — beginning of compression; 1=2 
— moment of formation of g black hole; 1=3 — end of 
calculation. 
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Figure 1. The occurrence of a shock wave. The Lagrange radius 
R Is applied along the axis of the abclssa, the 
logarlthn of density along the ordinate axis. Dif- 
ferent curves correspond to different moments of time. 
Here R_ = 1, A = 0.1, 

[Commas In this figure and In the following figures are equivalent 
to decimal points.] 









Oocurrehae of. i sonic wave, coming from the center 
= 0 • 75 j A = 0.5*R . PBH do not form. 






Figure 3. The same as in Figure 2 for R 
PBH do not fo’ n. 








Figure 5. Criteria of PBH. The Lagrange radius is applied on 
the axis of thp ahclssa and the ratio m/r on the 
ordinate axis. The moment of formation of a PBH 
corresponds to the tangency of the curve with a 
dashed horizontal. 











Figure 9. Conditions for the occurrence of PBH. R is applied 

• niu.x 

along the axis of the abcissa and A/R on the ordinate 
axis. All of the calculated variations are presented in 
the drawing., "0" — means the PBH does not form, " + " — 
PBH do form. The dashed curve shows the approximate 
boundary of the field of parameters during which PBH 
occur. 
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